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Higher orders in perturbation theory require the calculation of Feynman integrals at multiple loops. We 
report on an approach to systematically solve Feynman integrals by means of symbolic summation and discuss 
the underlying algorithms. Examples such as the non-planar vertex at two loops, or integrals from the recent 
calculation of the three-loop QCD corrections to structure functions in deep-inelastic scattering are given. 


1. INTRODUCTION 

Symbolic summation amounts to finding a 
closed-form expression for a given sum or se¬ 
ries. Systematic studies have been pioneered by 
Euler PP, and for specific sums, exact formulae 
have been known for a long time. Today, general 
classes of sums, for example harmonic sums, have 
been investigated (see e.g. Refs. [2]) and symbolic 
summation has further advanced through the de¬ 
velopment of algorithms suitable for computer al¬ 
gebra systems. Here, the possibility to obtain ex¬ 
act solutions by means of recursive methods has 
lead to significant progress, for instance in the 
summation of rational or hypergeometric series, 
see e.g. Ref. 0. 

In quantum field theory, higher-order correc¬ 
tions in perturbation theory require the evalua¬ 
tion of Feynman diagrams, which describe real 
and virtual particles in a given scattering pro¬ 
cess. In mathematical terms, Feynman diagrams 
are given as integrals over the loop momenta of 
the associated particle propagators. These in¬ 
tegrals may depend on multiple scales and are 
usually divergent, thus requiring some regulariza¬ 
tion. The standard choice is dimensional regular¬ 
ization, i.e. an analytical continuation of the di¬ 
mensions of space-time from 4 to D, which keeps 
underlying gauge symmetries manifest. Analyt¬ 
ical expressions for Feynman integrals in D di- 
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mensions may lead to transcendental or gener¬ 
alized hypergeometric functions, which have a 
series representation through nested sums with 
symbolic arguments. The main computational 
task is then to obtain the Laurent series upon 
expansion of the relevant functions in the small 
parameter e = (D — 4)/2. 


2. ALGORITHMS 


The basic recursive definition of nested sums is 
given by ^ 

S{n]mi, ...,mk]Xi, ...,Xk) = 
n 


where generally all \xi\ < 1. The sum of all 
rrii is called the weight of the sum, while the in¬ 
dex k denotes the depth. This definition actually 
includes as special cases the classical polyloga¬ 
rithms, Nielsen functions, multiple and harmonic 
polylogarithms piHTz] in their series representa¬ 
tions. For all Xi = 1, the above definition re¬ 
duces to harmonic sums |1I8I!II1()| and, if addi¬ 
tionally the upper summation boundary n —> oo, 
one recovers the (multiple) zeta values associated 
to Riemann’s zeta-function |2]. 

As an important property the 5'-sums in Eq. o 
obey the well-known algebra of multiplication. 
Specifically, any product 


S'(n;mi,...,m^;xi,...,a:*) 

xS{n;m[, ...,m'i;x[, ...,x'i), 


( 2 ) 




can be expressed again as a sum of single nested 
sums, hence in a canonical form, which is an im¬ 
portant feature for practical applications. The 
underlying algebraic structure in Eq. @ is a Hopf 
algebra, being realized as a quasi-shuffle algebra 
here, see e.g. Refs. [ 4111112111!] . The algorithm 
can be implemented very efficiently on a com¬ 
puter, see e.g. Refs. uninj. 

For the manipulation of the S-sums, we clas¬ 
sify certain types of transcendental sums. All 
sums in these classes can be solved recursively, 
i.e. they can be expressed in canonical form. The 
underlying algorithms realize a creative telescop¬ 
ing. They either reduce successively the depth or 
the weight of the inner sum, so that eventually 
the inner nestings vanish. Finally, the results can 
be written in the basis of Eq. o, ( as S'-sums with 
upper summation limit n) or as multiple polylog¬ 
arithms (which are S'-sums to infinity). 

Besides the quasi-shuffle algebra of multiplica¬ 
tion in Eq. the procedure relies on algebraic 
manipulations, such as partial fractioning of de¬ 
nominators, shifts of the summation ranges and 
synchronization of summation boundaries of the 
individual sums. 

Specifically, we consider convolutions. 
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(x' 

^7— ...,m'i]x' 2 , .■.,x'i ), 
conjugations. 


In all cases, the upper summation boundary 
should be consistent with the defining range of 
the binomials and the S'-sums. 

3. APPLICATIONS 

In perturbation theory, one way to classify 
Feynman integrals is according to the number of 
scales, i.e. the number of non-vanishing scalar 
products of external momenta or particle masses. 
According to this criterion, analytical expressions 
in D dimensions for the Laurent series in e either 
lead to transcendental numbers like (multiple) 
zeta values or (multiple) poly logarithms. Other 
classification criteria are, of course, the topology 
of a Feynman integral (number of loops and ex¬ 
ternal legs). 

3.1. One-scale problems 

Prominent examples of one-scale problems are 
massless two-point functions ll5l1bllVI . In par¬ 
ticular, the massless two-loop self-energy T1 has 
not only been of practical importance from a 
phenomenological perspective, but received also 
quite some interest from number theorists. For 
arbitrary powers of propgators, it is given by 

= y (Spy (Syj 

where = p 2 - q, p^ = pi - g, p^ = p^ - p^. 
Graphically, it is displayed in Fig.^ 
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Figure 1. The two-loop self-energy Tl. 

Here the interesting question has been which 
types of (transcendental) numbers appear in the 
e-expansion of this integral. For powers of the 
propagators of the form Ui = 1-1-e, i = 1,...,5 
it was known from explicit calculations up to the 
e®-term (heavily relying on symmetry properties) 














that multiple zeta values occur m However, 
it was unclear whether this suffices to all orders 
in e. Eventually, by deriving a double sum of 
(generalized) hypergeometric type, it was proven 
that multiple zeta values are indeed sufficient m 

Currently, with the help of symbolic summa¬ 
tion the e-expansion is known to the e^^-term CHI. 
The depth is limited by the fact that harmonic 
sums in infinity are expressed in a basis of tran¬ 
scendental numbers only up to weight 16. 

Another example for a one-scale problem, that 
received attention recently fTTinni is the non- 
planar vertex at two loops, which enters in calcu¬ 
lations of the quark and gluon form factors 
in QCD. Here the basic integral (displayed in 
Fig.m is given by 

X. . X /■ f d^P2 

where = p4 - p 7 , p^ = p^ - g, p^ = p 2 - pi, 
P 8 =P 2 - P3 and q'^ = {ps - p^Y. 



Figure 2. The non-planar two-loop vertex Vno- 

For general powers of propagators, Vno can be 
written as a double sum over Gamma functions, 
and if all = 1 , an expression in terms of hyper¬ 
geometric functions 3 F 2 and 4 F 3 has been given in 
Ref. Cni . After expansion, the sum can be solved 
in terms of the Riemann zeta function to any or¬ 
der in e using the algorithms for harmonic sums 
cm coded, as all our symbolic manipulations, 
in m We find the following expansion to order 

V^No(l, 1 , 1 , 1 , 1 , 1 ) = 

c 2 ( 2 X“ 2 - 2 e 

St (-9 ) 


1 5C2 , 27C3 

p4 ,2 , 


+23C2^ — e(48C2C3 — ll^Cs) + £^(-^( 2 ^ 
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-f I 2 I 8 C 2 C 7 + 1626C3 


20777 
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where we have taken out the usual MS-scheme 
factor 


( 477 )“^+'" 
" r(i-e) ■ 


( 9 ) 


This completes the section of examples with 
one-scale. 


3.2. Two-scale problems 

Nice examples of two-scale problems are pro¬ 
vided by the recent calculation of the three loop 
corrections in Quantum Chromodynamics (QCD) 
to the structure functions of deep-inelastic scat¬ 
tering 12.312412,^1 . Here, the two scales are the 
virtuality of the exchanged gauge boson = 
—q^ and the scalar product of the boson’s and 
nucleon’s momenta, 2 p ■ q, both combining to 
Bjorken’s dimensionless variable x = l{2p ■ q). 

The Feynman integrals under consideration can 
be expressed in nested sums and solved with the 
help of symbolic summation as follows. Imagine 
a mapping of a given integral I (x) depending on 
X, 0 < X < 1 to the space of discrete variables 
I{N), N G N, which is accomplished by means 
of an integral transformation, e.g. a Mellin trans¬ 
formation. Then one can obtain difference equa¬ 
tions for the Feynman integral 7(fV), which may 
be written as m 

ao{N) liN) + aYN) I{N - 1) + ... 

+amiN)I{N-m) = GiN), (10) 

where G{N) is some inhomogeneous term and Oi 
are some coefficients depending on N (and per¬ 
haps on e). The solution of Eq. cni) needs m 
boundary conditions 1 ( 0 ),..., /(m — 1). 

Single-step difference equations can be summed 
up analytically in closed form. Suppose we have 



















the equation 


ao{N) I{N) - ai{N) I{N - 1) = G{N), (11) 


then its solution will be 




n^iao(j) 
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( 12 ) 


In the case that the functions can be factorized 
in linear polynomials of the type N + m + ne 
with m, n being integer and N being symbolic, 
the products can be written as combinations of 
Gamma functions. In the presence of parametric 
dependence on e the Gamma functions should be 
expanded around e = 0, leading to factorials and 
harmonic sums. If the function G{N) is expressed 
as a Laurent series in e with the coefficients being 
combinations of harmonic sums in TV + m and 
powers of + m, m being a fixed integer, the 
sum in Eq. 03 can be done and I{N) will be 
a combination of harmonic sums in N + k and 
powers oi N + k with k being a fixed integer. 

Eq. n OH is an example of a recursion for Feyn¬ 
man integrals with dependence on symbolic pa¬ 
rameters. Solutions such as Eq. allow for 
an efficient implementation in computer algebra 
systems like m resulting in a largely auto¬ 
matic build-up of nested sums. For calculation 
of QGD corrections to structure functions men¬ 
tioned above, a systematic evaluation of nested 
sums was required for all integrals occurring in 
approximately 10.000 Feynman diagrams. Be¬ 
cause of the expressions being of excessive size at 
intermediate stages this task was well suited for 
the computer algebra system and the Summer 
package m for nested sums. 


3.3. Multi-scale problems 

Multi-scale problems arise in the calculation 
of cross sections with more kinematical invari¬ 
ants, like e.g. jet cross sections. The methods 
and algorithms for generalized sums have already 
been used in full-fledged QGD calculations, for 
instance in the evaluation of higher order correc¬ 
tions to e+e“ —> 3jets m 


In general, Eqs. may also be used to 

expand higher transcendental functions in a small 
parameter around integer values. Starting from 
the series representation of, e.g. the first Appell 
function 


Fi{a,bi,b2-,c;xi,X2) = 


OO OO 
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I^mi+m 2 rriil TO 2 ! ’ 
or the second Appell function 
F2{a,bi,b2;ci,C2;xi,X2) = 


OO 

E 

7711=0 7712=0 


E 




(13) 


(14) 


we see that Eqs. (0-© apply if the expansion 
parameter e occurs in the argument of the rising 
factorials (Pochhammer symbols), defined as 


= f(^x)f{x+l)...f{x + m-l). (15) 


It should also be stressed at this point, that al¬ 
though the definition of the S'-sums in Eq. o is 
very general, the specific algorithms for convolu¬ 
tion, conjugation etc. are subject more restrictive 
assumptions. 

In particular, the algorithms underlying the 
evaluation of Eqs. do rely on the fact 

that the modulus of the summation index in 
the argument of the S-sums or in the denomina¬ 
tors is always one. This changes, if for instance 
hypergeometric functions jFj-i (or more gen¬ 
erally Gamma-functions) are expanded around 
half-integer values. Such a situation occurs for 
example in the calculation of massive higher loop 
integrals in Bhabha scattering |23|. 

Some extensions of the summation algorithms 
for expansion around rational values, leading e.g. 
to binomial sums and inverse binomial sums are 
discussed in Ref. |29I30| . 


4. CONCLUSION 

Symbolic summation has advanced to an im¬ 
portant method for the calculation of higher order 
corrections in perturbative quantum field theory. 
The field has seen significant progress during the 
past years and we have given various examples 










from complete calculations, e.g. the recent eval¬ 
uation of the third-order contributions in pertur¬ 
bative QCD to the structure functions of deep- 
inelastic scattering ESE3ES1- These cutting edge 
calculations show that the method of symbolic 
summation provides very powerful means for the 
practical computations of Feynman diagrams. 

In closing, we note that all practical applica¬ 
tions do heavily rely on computer algebra imple¬ 
mentations of the algorithms discussed here. In 
the symbolic manipulation program m, which 
is a fast and efficient computer algebra system 
to handle large expressions, harmonic sums can 
be manipulated with the Summer package cni- 
For the S'-sums of Eq. m there exists an exten¬ 
sion, the XSuMMER package m in , which im¬ 
plements algorithms of Eqs. 0-©. As an al¬ 
ternative within the GiNaC framework m the 
package nestedsums m provides similar func¬ 
tionalities. Very recently, also Ref. m appeared, 
which limits itself to the problem of expanding 
hypergeometric functions jFj-i around integer 
parameters to arbitrary order and provides an im¬ 
plementation in Mathematica . 

We believe all these packages may also be useful 
for a larger community. 
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